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ABSTRACT 

We present a study of dense structures in the L 1495 filament in the Taurus Molecular Cloud 
and examine its star-forming properties. In particular we construct a dust extinction map of the 
filament using deep near-infrared observations, exposing its small-scale structure in unprecedented 
detail. The filament shows highly fragmented substructures and a high mass-per-length value of 
■Miine = 17 Mq pc -1 , reflecting star-forming potential in all parts of it. However, a part of the filament, 
namely B211, is remarkably devoid of young stellar objects. We argue that in this region the initial 
filament collapse and fragmentation is still taking place and star formation is yet to occur. In the 
star-forming part of the filament, we identify 39 cores with masses from 0.4 . . . 10 M Q and preferred 
separations in agreement with the local Jeans length. Most of these cores exceed the Bonnor-Ebert 
critical mass, and are therefore likely to collapse and form stars. The Dense Core Mass Function 
follows a power law with exponent V = 1.2 ± 0.2, a form commonly observed in star-forming regions. 
Subject headings: ISM: clouds — (ISM:) dust, extinction — ISM: individual (L1495) — ISM: structure 
— stars: formation 



1. INTRODUCTION 

Filaments appear to be common structural features 
in both quiescent and star-forming molecular clouds. 
In particular, r ecent observations of molecular clouds 
with Herschel (jPilbratt et al.l 120101 ) have signified the 
role of filaments as a momentous, perhaps even domi- 
nant mode of star f o rmation (e.g. iHenning et al.l 120101 : 
iMolinari et "all 120101: iMen'shchikov et al.l 120101 ) . How- 
ever, their formation and small scale struct ure is still not 
well understood and und er debate (see e.g. IMversI 120091 
and references therein). ISchneider fc Elmegreenl (|1979f ) 
noted that filaments tend to fragment into equally spaced 
condensations of sizes ~ 10° pc. These clumps, in turn, 
consist of even smaller, gravitationally bound entities, 
namely cores, of sizes ~ 10 _1 pc, which is of the order of 
the local Jeans length. These cores are then believed to 
be the dir ect precursors of stars a nd binary systems (see 
review bv lBergin fc Tafallall2007l) . 

Such a hierarchical structure can also be found in the 
Taurus Molecular Cloud (Fig.[T|), which is one of the clos- 
est star-forming reg ions at a distance of only 137 ± 10 pc 
(|Torres et al.l[200 7), and therefore an excellent testbed 
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to study the small sc ale structure of filaments. Large 
scale surveys in CO (Mizuno et all 1995; O nishi et al.l 
[1991 INaravanan et aLll2008D and mid-infrared dust emis- 
sion dRebull et al.l 120101), alongside with dust extinc- 
tion maps ( Dobashi ct al. 2005; Kainulainen et al.| [2009b; 



lLombardi et al.l I2010D have built a detailed view of 
a complex network of filaments in Taurus, which are 
known to be the birth sites of st ars (see reviews about 
star formation in T aurus by e.g., IPalla fc Stahlerl 120021 : 
iKenvon et al.l I2008D . These filament are composed o f 
clumps (|Barnardlll927tlLvndsl[T962l:lLee fc Mverslll999D 
with sizes of ~ 1 pc, which host a plethora of cores with 
typical sizes of ~ 0.1 pc. These have been detected e.g. 
via molecular line observations tracing cold and dense 
gas with volume densities of n H2 > 10 5 cm~ 3 , such 
as N H 3 (IJiiina et al.l Il999l) and H 13 CO+ (|Onishi et al.l 
2002, hereafter |OMKQ2). Further down in this hier- 
archical structure, from dust continuum observations 
ISadavov et al.l (|2010f > found a population of pre- and pro- 
tostellar cores, supposedly representing the scale of direct 
core-protostar connection. The Taurus Molecular Cloud 
is known to host > 250 Young Stellar Objects (YSOs, 
iRebull et al.ll20ldh . and their distribution shows a strong 
correlation to the filamentary structures, i.e. their birth 
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Fig. 1. — Extinction map of the Taurus Molecular Cloud (Kainulaincn ct al. 20093) deri ved from 2M ASS data. Contours are plotted 
in steps of Ay = 2 mag. The circle marks the center of the hub-filament system in L 1495 (Myers 2009) with three filaments emanating 
towards the east, south-west and west. The thick solid lines mark the region of our Omega2000 observations in L 1495, and the control 
field to the south of it. 



sites (|Hartmannl 120021: iKenvon et al.ll2008l) . 

To further understand the connection between fila- 
mentary structures and star formation, it is essential to 
characterize the filament in high spatial resolution over 
a wide dynamic range. Molecular line species can be 
used as probes for column densities, but typically rather 
exhibit a narrow dynamical range. Larger dynamical 
range can be achieved by thermal dust emission obser- 
vations, which are used to trace intermediate to high 
column densities, and dust extinction mapping, which 
is more sensitive a t low to intermediate column den- 
sity range (see e.g. iPineda et al.ll2008t iGoodman et al.l 
2009; iVasvuninaet al.l 120091 for a comparison of differ- 
ent methods). The aforementioned observations allow 
to determine the structure of the low density clumps on 
one hand, and the core population detected by dense gas 
tracers on the other hand. However, to make a direct 
connection between the diffuse envelope structure and 
the denser cores, one preferably needs a single tracer with 
a uniform calibration, which covers both these regimes. 

Therefore, the aim of thi s work is to create a col- 
umn density map of L 1495 (|Lvndslll962l ) with unprece- 
dented dynamical range and resolution with the aid of 
near-infrared (NIR) dust extinction mapping. NIR col- 
ors of stars in the background of clouds provide hun- 
dreds to thousands of pencil beam measurements of red- 
dening, which can be smoothed out to construct maps 
of column density along a molecular cloud. This is the 
basis of NIR excess technique nice and its derivatives 
(jLada et al.lll994t ILombardi fc Alvesl 120011: Hbster et al.l 



120081 ILom bardi 2009) , which obtain extinction measure- 
ments by correlating the observations from the science 
field with a nearby control field. The color excess with 
respect to this control field is then completely attributed 
to be caused by the molecular cloud. These methods 
typically probe extinctions Ay ~ 0. 5.. .50 ma g (e.g . 
I Lombardi fc Alvesl 120011: iKainulainen et al.l [20061 120071; 
iRoman-Zuniga et al.ll2009l ). allowing us to examine both 
the detailed fragmentation of the filament into dense 
cores, and the envelopes surrounding them. Thus, the 
observations presented in this paper will provide the most 
complete census of the small scale structures in the Tau- 
rus filament so far, and shed more light on the link be- 
tween filaments and star-forming cores. 

The paper is organized as follows. In we give an 
overview of the observations and data reduction. In 
we present our NIR extinction map of unprecedented 
depth and resolution, and derive the dense core popu- 
lation, which is discussed in section SjJJ Finally, the con- 
clusions and summary follow in Sj5l 

2. OBSERVATIONS AND DATA REDUCTION 

We carried out NIR observations of the Taurus fil- 
ament L 1495 with the Omega2000 camera at the 
Calar Alto 3.5m telescope. The camera is equipped with 
a HAWAII2 HgCdTe detector that offers 2048x2048 pix- 
els with a pixel scale of 0.45 " and a field of view of 
15.4 x 15.4'. Observations were carried out during two 
observing periods f rom October 2004 to December 2005 
(|Quanz et al J 12010) and January to March 2009. In to- 
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Fig. 2. — Our star-galaxy separation is based on the stcllarity 
index C, elongation E and its normalized FWHM Fq. The hatched 
areas indicate the regions populated by stars. Due to illustrative 
reasons, only a small part of our final catalog is shown. 
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Fig. 3. — Photometric uncertainties for our sources in the differ- 
ent filters. They contain contributions from the photometric er- 
ror and the uncertainty of their individual photometric zero-point. 
This figure clearly shows the contributions from different frames, 
which differ in terms of sensitivity due to varying observing con- 
ditions. Due to illustrative reasons, only a small part of our final 
catalog is shown. 



tal, 21 science fields and 1 control field were observed in 
J, H, and K s bands (1.2 //m, 1.6 /^m, and 2.2 /im, respec- 
tively). The locations of the observed frames are shown 
in Fig.[2 

The observations of a single Omega2000 field consisted 



of 30 dithered exposures, each having an integration time 
of 60 s. This resulted in total exposure times of 30 min- 
utes per filter per field. We performed a st andard NIR 
data reduction with the Omega2000 pipeline (Fass benderl 
|2003[ K including Hatfield, dark, and bad pixel corrections, 
sky subtraction, and image registering. The world coor- 
dinate s ystem (WCS) was appended to the frames with 
koords (|Goochl I1996D by matching the positions of at 
least 10 stars per frame to matching stars from the cor- 
responding 2MASS image. This resulted to the rms ac- 
curacy of 0.1 px (or 0.045 ") in all frames. 

In addition to stars, a significant amount of back- 
ground galaxies was expected to be detected in the 
frames. In order to disentangle these two types of sources 
from each other, we performed pho tometry using S Ex- 
tractor (jBertin &; Arnoutsl Il996[ h which provides an 
automatic source classification via the so-called stellarity 
index (CLASS_STAR). In particular, photometry was run 
in a fully automated two-pass mode. First, SExtrac- 
TOR was run for the brightest stars only, in order to get 
an accurate estimate of the seeing, which is a crucial in- 
put parameter to get reliable measures of the stellarity 
index, and furthermore determined our choice of aper- 
ture size. Second, the actual photometry was performed, 
resulting in instrumental magnitude and uncertainty, po- 
sition, the stellarity index C, elongation E (the inverse 
of the aspect ratio) and the parameter F$, which is a 
source's FWHM normalized by the seeing. 

Each source was classified as either star or galaxy 
follow ing the classification scheme from iCantiello et al.1 
(2005) . In this scheme, a source is regarded as a star if 
it fulfills two out of the three following conditions: 

• C > 0.5 

• E > 0.8 

• 0.75 < F < 1.25 

which are indicated as hatched regions in in Fig. [2] 

The photometric calibration was done for each frame 
separately by matching sources within the frame with 
stars from the 2MASS catalog, and determining photo- 
metric zero-points and color corrections using those stars. 
These catalogs in J, H and K s were then merged, and 
sources with a detection in only one filter were rejected. 
Saturated stars were replaced by their 2MASS coun- 
terparts. The final catalog of L 1495 contains ~33,000 
sources, of which 63% are classified as stars, according 
to our classification scheme. Photometric uncertainties 
for an excerpt from this catalog are shown in Fig-El 
The median magnitudes for which a = 0.1 mag are at 
J = 20.6 mag, H = 19.4 mag and K s — 18.8 mag, which 
is on average 4 mag deeper than 2MASS, and emphasizes 
the depth of our observations. 

3. RESULTS 
3.1. Dust Extinction in the Taurus Filament 

We used the JHK S photometry of the detected sources 
to derive dust extinction through the Taurus filament 
L 1495. In the following, we first shortly describe the 
adopted technique and then present the results. 

As mentioned earlier, in addition to stars there is a sig- 
nificant number of background galaxies among detected 
sources. The intrinsic colors of galaxies are very different 
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Fig. 4. — NIR color-color diagram of the science (left) and the 
control field (right). A shift cause d by an extinction of Ay = 5 mag 
(assuming the reddening law of Rieke & Lcbofsky 1985) is indi- 
cated by the arrow in the left panel. T he continuous line shows an 
unrcddcncd Zero-Age Main-Sequence (Sicss et al. 2000) for com- 
parison. 



from those of stars, and therefore theses two classes need 
to be treated separately in deriving extinction towards 
them. This difference is illustrated in Fig.[H which shows 
color-color diagrams of our reference (right panel) and 
one science field (left panel). Especially in the reference 
field, one can clearly see the different color distribution of 
stars and galaxies, as it is unaffected by the dust extinc- 
tion from the molecular cloud. We note that there is a 
concentration of stars at [H — K s , J — H] = [0.7, 1.0] that 
is clearly outside the Zero-Age Main-Sequence (ZAMS) 
of stars. These sources could be YSOs, brown dwarfs, or 
even point l ike distant galaxies and or quasars. 

Recently, iFoster et"all (|2008| ) presented a NIR color- 
excess mapping technique tailored specifically for inclu- 
sion of galaxies. In this technique, the colors of stars 
behind the cloud are compared to the intrinsic colors of 
stars in the reference field. In addition, the magnitude 
dependent colors of galaxies are compared to those in the 
reference field. This comparison yields a pencil-beam-likc 
extinction value towards eac h source (for further details 
we refer to lFoster et al.ll2008[ ). From these values we then 
generate a uniformly gridded extinction map by smooth- 
ing the data using a Gaussian kernel. Sigma-clipping re- 
moves sources with outliers in extinction, which might 
be caused by intrinsically different colors (e.g. YSOs, 
brown dwarfs, etc). For our map, we chose a resolu- 
tion of 0.9 ' as a compromise between resolving the small 
scale structures in the filament and having low noise 
over the dynamical range covering most of the map. In 
this resolution, the map covered the dynamical range of 
Ay = 1.5 ... 35 mag. The noise of the final map depends 
on extinction, being a = 0.5 mag at Ay = Omag and 
cr = 2 mag at Ay = 35 mag. 

Fig. [5] shows the resulting extinction map, with blow- 
ups of different regions of it shown in Fig. [6] The map 
reveals, on one hand, a remarkably well defined, high 
aspect ratio filamentary structure, and on the other 
hand complex clumpy substructures within this filament. 
We note that within the mapped region we find six 
Barnard objects, namely the clumps B 211, B 213, B 216, 
B 217, and B 218 (which represent the filament), and B 10 



TABLE 1 

Barnard Objects in the L 1495 filament 



Region 


Mass 


Area 






[M ] 


[pc 2 ] 


[M pc" 2 ] 


B211 


84 ±5 


0.40 ± 0.06 


210.7 ± 14.8 


B213 


87 ±6 


0.48 ±0.07 


180.3 ± 12.6 


B216 


57 ±4 


0.34 ±0.05 


169.8 ± 11.9 


B217 


29 ±2 


0.19 ±0.03 


154.8 ± 10.8 


B218 


22 ± 1 


0.13 ±0.02 


174.8 ± 12.2 


Total 


280 ± 19 


1.54 ±0.22 


182.3 ± 12.8 



Note. — All parameters are calculated for re- 
gions with Ay > 5 mag. 

a Values for mass and area in the table are 
rounded, and therefore their division leads to dif- 
ferent values compared to the correct calculation. 

(which forms the hub). 

Within the filament, the lowest completely mapped iso- 
contour is at Ay — 5 mag. The total mass within this 
region is 280 ± 19 M Q , yielding a mass surface density of 
Sm = 182 ± 13M Q pc -2 . TableQ] lists these parameters 
also separately for the subregions. The mass was esti- 
mated from the extinction map by assuming all hydrogen 
in molecular form and the standard conversion factor of 
NvJAy = 0.94 x 10 21 cm^rnag- 1 (IBohlin et al.lll978l: 
iRieke fc "T cbofskv 1985), which then yields the total hy- 



drogen mass 



M 



= 4.89 x 10" 



D 



137pc 



0.45' 



.4 



1 mag 



(1) 

where D is the distance to the cloud, 9 is the pixel scale 
of the map, and Ay^ are the extinction values of the 
individual pixels. Using a standard clou d composition of 
63% hydrogen, 36% helium and 1% dust (jLombardi et al.1 
[2001 one can assume the total mass from Equation ([TJ 
to be 

M tot = 1.37M H . (2) 

The blow-ups shown in Fig. [5] reveal numerous dense 
fragments with peak extinction values of Ay > 15 mag. 
This population of dense core£\ will be investigated in 
more detail in the following section. The main filament 
is strongly meandering along almost its entire length of 
8 pc, and sometimes forms even ring-like structures (e.g. 
B216, B10). However, B211 is different in that it ex- 
hibits a linear, narrow bar with a length of 1.5 pc and 
a median FWHM of only 0.11 pc. The inter-core re- 
gions also show considerably higher extinction values of 
^ 10 mag in contrast to the other parts of the fila- 
ment, where cores seem to be rather isolated. 

3.2. The Dense Core Population 

Identification and characterization of small-scale struc- 
tures, such as cores, from 2D data (e.g. extinction 
maps or continuum observations) is a long-lasting, 
non-trivial problem, and several methods exis t (e.g., 
IStutzki fc Guested Il990l IWilliams et all [l99l . We 

1 Our analysis ( ^4.21 1 shows that these are indeed dense cores, i.e. 
gravitationally bound entities with densities of njj 2 ~ 10 4 cm - 3 . 
Therefore, we already introduce this term here, for reasons of con- 
sistency. 
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Fig. 5. — Extinction map of the Filament L 1495 with a resolution of 0.9' derived from deep NIR observations with Omega2000. 
Contours are plotted in steps of A\r = 5 mag. We separate the filament into different subregions, which conform to Barnard's Dark Objects 
(Barnard ct al. 1927). The boxes indicate the positions of the zoom-ins shown in Fig.[6] 



adopted a two-step approach, which consisted of back- 
ground removal and 2D thresholding. In the first step, 
we performed a subtraction of large-scale structures from 
the map using wavelet filtering. This resulted in a map 
describing structures only at s patial scales smal ler than 
0.14pc ("cores-only" map, c.f. lAlves et al.ll2007f ). In the 
second step, we used the 2-dimensional thresh olding al- 
gorithm Clumpfind2D (jWilliams et all IT994T ) to iden- 
tify a population of dense cores from the cores-only map. 
The algorithm works by first contouring the data (usu- 
ally these levels are defined by a detection threshold, 
and subsequent contours with separations of multiples 
of the rms noise of the observations) , and then searches 
for peaks of emission, which mark the cores. It then 
follows these down to lower intensities, assigning every 
pixel to a certain core (if any). Obv iously, such an algo - 
rithm works best for isolated cores. iPineda et all (2009) 
pointed out that for observations, where one has to deal 
with blending of cores, CLUMPFIND2D is very sensitive 
to the choice of input parameters. Therefore, we ran 
a set of calculations with different values of detection 
thresholds and level separations. We performed various 
tests, which included visual verification of the positions, 
sizes and shapes of the cores. A detection threshold and 
a level separation of 3<r was chosen to yield a represen- 
tative core dis tribution. This is i n agre ement with the 
simulations by iKainulainen et all (|2009al ) , in which such 
parameters resulted in an acceptable core identification. 

However, the narrow bar in B211 was not removed 
by the wavelet filter, eventually not allowing a look at 
the "naked cores". This prohibited a reliable recovery 
of the dense cores and their properties in B211. We be- 
lieve that this is caused by its different structure, which 
is discussed in more detail in § 14.11 Therefore, we did 
not consider this region in the following investigation. 
Furthermore, we also did not take into account B 10 be- 
cause of its cluster-like nature and very high extinction, 
that does not allow to break the st ructure into sepa- 
rate cores (jKainulainen et ail l2009afl . Follow-u p anal- 
ysis using dendrograms |Rosolowsk*v et~afll20"08h on 3D 
position-position-velocity molecular line data might al- 
low us to also disentangle the dense core population in 



these two regions in the western part (L 1495-W). In this 
paper, however, due to the limitations of 2D data we re- 
strict our analysis to the dense core population in the 
eastern part (L 1495-E), which consists of regions B 213, 
B216, B217, and B 218. 

In total 39 cores were detected in L 1495-E (Fig. Ta- 
ble[2]). We calculated the core orientation together with 
the zeroth, first and second moment, which after suit- 
able conversions conformed to the core mass, position 
and FWHM along the minor and major axis. Further- 
more, we defined the effective radius as r c ft = (A/ir) 1 / 2 , 
assuming spherical symmetry. We note that the mass 
is quite insensitive for the choice of the threshold level, 
changing on average only 15% when going from 3er to 6a. 
A few cores (#23, #29 and #35) extended beyond our 
mapped region. We used data from the 2MASS catalog 
to extend the map to cover these cores completely. 

The nearest neighbor separation of the cores exhibits 
a strong peak at 0.13 pc with a dispersion of 0.04 pc, and 
all cores had at least one neighbor within 0.25 pc. An- 
other approach to quantify the "clusteri ng" is to measure 
the two-point correlat ion function fe.g.. [Johnstone et al.1 
120001 : lHartmannll2002[ ) . where the number of core pairs N 
with separations in the interval [log(r), log(r) + dlog(r)] 
is compared to the number Nq in a sample, where cores 
were randomly distributed over the whole region. The 
two-point correlation function is then 



( \ N 
W{r) = Ao 



1. 



Applying this to our sample of cores, strong correlations, 
i.e. w(r) > 0, on separations of r ~ 0.15 pc and r ~ 
5pc are evident (Fig. [7]). The latter, large separation 
derives from small clusters of dense cores at both ends of 
L 1495-E, whereas the first one represents the separation 
between the individual cores inside these clusters. At our 
pixel scale of the map, a separation of 0.15 pc corresponds 
to ~ 7px, and is, therefore, well above the resolution 
limit. 

Core masses range from M core = 0.4 . . . 10 M©. The 
total mass found in cores is 104 ± 16 M©, which makes 
up ^50% of the total mass in L 1495-E (as defined by 
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Fig. 6. — Zoom-ins into the extinction map. The scale bar indicat ed in the bottom p anel is com mon to a ll sub-plots. Contours are plotted 
in steps of 5o~ each . The positions of dense cores mapped in NH 3 Uiiina et al.lll999l '), H 13 CO+ (QMEM) are overplotted alongside YSOs 
HRebull et aUrSOlCf) . 



the Ay = 5 mag contour). We derived mean core den- 
sities p = 3M core /(4r 3 ff ), and hydrogen number den- 
sities nfj 2 = p/((imn), where mn is the proton mass 
and /i = 2.34 the mean molecular weight per hydro- 
gen molecule, assuming a molec ular cloud with standard 
composition (|Lada et al.l 120081 ). The frequency distri- 
bution of the derived number densities closely follows a 
Gaussian distribution with nn 2 = 14-3 x 10 3 cm~ 3 and a 
dispersion of 4.1 x 10 3 cm~ 3 . 

In Fig. [8j we show the mass-size relation of the cores 
in our sample. For reference, the figure also shows the 
sensitivity limit of Ay = 1.5 mag for an isolated, flat 
and spherical core. lLada et al"] (|2008l ) derived a slope 
of of M oc r^g for the dense cores of the Pipe Neb- 
ula. Due to the low number of cores and the strong 
influence of the completeness limit on the fit, we rather 
show representative relations with constant column den- 
sity [M oc f*^ ff ), whic h would be p redicted for clouds 
obeying Larson's laws (|Larsod fl981). and constant vol- 
ume density (M oc r 3 ff ), which would be expected in 
case of constant thermal pressure and kinetic tempera- 
ture. However, with the data at hand we cannot rule out 
one or the other model. 

The Dense Core Mass Function (DCMF), which is 
the number of cores per logarithmic mass interval 
diV/dlog(m) oc m~ r , exhibits its maximum at ~ 2.0 Mq, 
and falls off to both the low-mass and high-mass ends 
(Fig.O. This peak could be either true, or simply be 
caused due to incompleteness. However, estimating the 



completeness an d accuracy of the DCMF is not a triv- 
ial ta sk (see e.g., IKainulainen et al]|2009al IPineda et al.l 
2009). This is partly due to the inability of the adopted 
algorithm to detect cores from the variable, more ex- 
tended column density component. Partly it is also due 
to the noise in the map, which makes mass determina- 
tion of low-mass cores significantly less accurate than 
for high-mass cores. In an a nalysis similar to this pa- 
per, |Kainulainen^eF^l] (|2009aD analyzed the accuracy de- 
rived for the cores in the Pipe Nebula, and concluded the 
DCMF to be accurate above about 1.2 Mq. By compar- 
ing the noise levels within these two maps we roughly 
estimate the DCMF to be accurate for M core > 1.5 Mq. 
A power-law fit to the DCMF above this completeness 
limit yields an ex ponent T = 1.2 ±0.2. This slope agrees 
with the results of lOMK02t which found T = 1.5±0.3 for 
3 < M/M < 20. Slopes of T ~ 1.2 are als o commonly 
observed in other star-forming regions (e.g. iMotte et al.l 
[1991 [Rathborne et al.ll2009t lAndre et al.lf2oTo( )~ 

The distribution of aspect ratios R, which is defined 
as the ratio of the major vs. the minor axis a/b, is 
strongly skewed. The mode of the distribution is found 
to be R = 1.3, whereas the mean aspect ratio is found 
to be R = 1.5, with a dispersion of 0.4 (Fig.rTUl left 
panel). Within the region of our observations, we find 
five cores from the s ample of optically selected cores from 
iLee fc Mveri (fl999l ). Three of these cores were classi- 
fied as round (R, = 1), and the remaining two had only 
moderate aspect ratios of R ~ 1.8. Compared to the all- 
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Fig. 7. — The distribution of cores in L 1495— E. The core-IDs can be found in Table[2] In the top right corner, the inset depicts the 
two-point correlation function w(r). 




r eff [P c ] 

Fig. 8. — Mass-size relation for our cores. The hatched area 
marks the region of incompleteness. The dotted lines represent 
the cases of constant volume density M ex. r 3 and column density 
M oc r 2 . 



sky sample of 406 cores from iLee fc Mversl (|1999f ) with 
R = 2.4 ± 0.1, the cores in Taurus are therefore only 
moderately elongated. 

We also considered the orientation of the cores with 
respect to the filament. Due to its meandering nature, 
this is not a trivial task. Therefore, we first defined the 
baseline of the filament by setting points along its path, 
and then interpolated between these points using a cubic 
spline function. The orientation angle of the core with 
respect to this baseline was then obtained, and we find 
that the 39 cores seem to have a slight tendency to align 
with the filament. However, a KS-test results in a prob- 
ability of p = 0.46 that this distribution is drawn from a 
uniform distribution. Given the large error bars on the 
individual orientations, which come along with the small 
aspect ratios, we can therefore not draw any secure con- 




M [MJ 

Fig. 9.— The Dense Core Mass Function (DCMF) of L 1495- 
E, which was obtained by smoothing the core masses (indicated 
as vertical dashes along the abscissa) with a Gaussian kernel with 
FWHM=0.3dex. Points with error bars are plotted in separations 
of 50% of the FWHM. The dashed line shows the power law fit with 
slope T = 1.2 ± 0.2 for M > 2.OM . The hatched area indicates 
the region of incompleteness. 



elusions about preferential or random core alignment in 
this filament. 

3.3. Reddening Law 

Algorithms like GNICER obtain the excess extinction 
of a source by comparing its position in the color- 
color diagram with respect to a control field. It is 
crucial for a proper estimatio n of extinction values to 
assume a reddening la w (e.g., iRieke fc L cbofskv 119851 : 
llndebetouw et alj |2005). For observations in JHK S the 
ratio of the colour excesses 



a 



E(J - H) 
E(H - K s ) 



(3) 
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TABLE 2 
Dense Cores in L 1495-E 




1 1.5 2 2.5 -60 -30 30 60 

Aspect Ratio a/b Orientation [deg] 



Fig. 10. — left panel: Distribution of aspect ratios R, which was 
obtained by smoothing the data points (indicated as vertical dashes 
on the abscissa) with a Gaussian kernel of FWHM= 0.3. Data 
points are plotted in steps of 50% of the FWHM. — right panel: 
Distribution of core orientations with respect to the filament. The 
FWHM for the Gaussian kernel was chosen to be 30 ° . 

is constant, and defines the direction of the displacement 
of a source in the color-color diagram when it is subject 
to extinction. Therefore, the colors (J—H) and (H — K s ) 
are connected by the linear relation 

(J-H)=a(H- K s ) + (3. (4) 

As unreddened ZAMS stars occupy only a small region 
of the color-color-diagram, the slope a could be simply 
obtained by fitting Equation (U) for all stars simultane- 
ously. This, however, would give the high number of 
low-extinction stars too much weight compared to the 
low number of high-extincti on stars. Therefore, we fol- 
lowed a procedure similar to lLombardi et all (|2006t ). 

The first step was to obtain visual extinctions Ay of ev- 
ery star by converting color excesses E ( J—H) and E(H— 
K s ) a ssuming the reddening law of llndebetouw et all 
(200$. Stars are then binned in the color-color diagram 
according to their extinction Ay. These data points are 
then fitted according to Equation d4j to obtain a cor- 
rected slope. With the aid of this new slope, corrected 
extinction values Ay are then computed. This procedure 
is repeated until the slope a converges, which generally 
happened after 3-6 iterations. In order to decrease de- 
pendence on the bin size, we repeated this procedure 
for bin sizes ranging from A Ay = 0.5 mag ... 1.5 mag in 
steps of 0.2 mag. The best fit slope a was then the vari- 
ance weighted mean from all runs with different bin sizes. 

For this analysis we considered only sources classified 
as stars without any photometry flags. The best fit slope 
was found to be a = 1.82 ± 0.08. This value is in very 
good a greement to the reddeni ng laws of |Lombardi et all 
(2006, a = 1.82 ± 0.03) and llndebetouw et all (|2005l 
a = 1.78 ±0.15). 

4. DISCUSSION 

4.1. Star Formation in the Filament 

Due to its proximity, the Taurus Molecular Cloud is 
an excellent target to study the connection between fil- 
amentary structure and star formation. In this paper 
we present a new, high-resolution view of the L 1495 fil- 
ament, allowing characterization of its density structure 
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down to tenth-of-a-parsec scales. Indeed, as illustrated 
in Figs. [5][2 L 1495 shows a highly fragmented small- 
scale structure and in some places extinctions, which are 
higher than can be traced by the chosen technique (i.e. 
Ay > 35 mag) . 

To estimate the star-forming potential of the fila- 
ment, we calculated the mass-per- length Mi; no . For an 
undisturbed filament, filaments with a mass exceeding 
the critical value of Mj™* = 15M w pc~ * are prone to 
filamentary fragmentation and collapse ( Ostrikcr 119641 : 
llnutsuka fc Mivamal I1997D . In presence of magnetic 
fields, this critical value can be increased by a factor 
of 2 (|Fiege fc Pudritzll2000af) . Due to the presence of 
dense cores, L 1495 does not really fulfill the require- 
ment of being an undisturbed filament. However, con- 
sidering only the large scale structure, which was fil- 
tered out in the wavelet transform, we derived a value 
of Miino = 16.8 ± 0.2 Mq pc _1 for the whole filament. As 
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Fig. 11.— Probability Density Function (PDF) of B213 and 
B211. Our observations do not completely cover regions with 
Av < 5 mag (hatched area). The dashed line is the log- normal 
fit to B 213, which is shown as a scaled-down version in the right 
panel. 

this must be seen as a lower limit, as some of the mass is 
already bound in the dense cores, this value underlines 
the star- forming potential of the filament. 

In Fig. [6] we show the previously known YSOt0. Inter- 
estingly, while YSOs are generally found close to dense 
cores in the fil ament, the bar-like B 21 1 is completely de- 
void of them ((Goldsmi th et al.l l2008f ) . Given that the 
mass-per-length value in this particular region is also 
high (Afiino = ISMqpc^ 1 ), and that it also harbors sev- 
eral H 13 CO + cores signifying the presence of high density 
gas (n cr j t = 10 5 cm~ 3 ), it can be argued to be well un- 
derway towards star formation. 

More insight into the assembly of material in the non- 
star-forming part of the filament can be gained by ex- 
amining the distribution of column densities. Therefore, 
in Fig.QT] we compare the quiescent B211 with B213, 
in which star formation has already occurred. The col- 
umn density distribution of B213 shows a continuous 
decline towards higher extinctions. In contrast, B211 
shows a flat plateau at intermediate extinctions, falling 
steeply at extinctions higher than Ay ~ 13 mag. This 
excess with respect to the column density distribution of 
B213 can be related particularly to inter-core regions. 
This is illustrated in Fig.fl^J which shows the contour of 
Ay = 11 mag for both B 211 and B 213. Clearly, in B 213 
such extinction range is related to cores, while in B 211 it 
is more continuous and related to the material surround- 
ing the cores. Therefore, we suggest that B211 is not 
at all devoid of star-forming capability, but in an early 
stage of its evolution where the filament is fragmenting 
and the inter-core regions have not yet been cleared from 
material by accretion towards the cores. Given the high 
mass reservoir in the filament, such fragmentation is in- 
evitably leading to star formation in B 211. 

In L 1495-E, the star-forming part of the filament, we 
investigated the spatial distribution of YSOs with re- 
spect to the filament. Class I and Flat-Spectrum sources 
are found at extinctions of Ay = 8.2 ± 3.0 mag, which 
is higher compared to the extinctions at the positions 

2 IRebull et al. 1 12010'') separate YSOs into Class I, Flat-Spectrum, 
Class II and Class III. We include all YSOs, which were classified as 
confirmed, probable, and possible Taurus members in their paper. 



B211 




Fig. 12. — Zoom-ins into the regions B211 and B 213 with con- 
tours at Ay = 5 mag (gray) and Ay = 11 mag (black) . 



1 


G 
O 




% • 




o 

• • 


non— equilibrium 


8««..o. 

• o • • 


. equilibrium 

• 


• • 
•• 


• 


• without 

OMK02 counterpart 

O with 


• i 





0.5 1 2 5 10 



Fig. 13. — Ratio of core masses Al corc to the Bonnor-Ebert crit- 
ical mass Mbb plotted against core mass. Mbe was calculated 
individually for every core, assuming its respective density and a 
temperature of T = 10 K. Open circles repre sent core s which have 
a counterpart from the dense core sample of O MK02t 

of Class II and III (Ay = 5.2 ± 2.6 mag). The mean 
spatial separation of Class I and Flat-Spectrum sources 
from the filament was found to be only 0.06 pc, whereas 
older sources showed considerably larger mean distances. 
To verify whether this could be a pure evolutionary ef- 
fect, we performed a simple Monte Carlo simulation to 
model the time dependence of the spatial distribution of 
sources. We distributed 10,000 test sources with a 3D ve- 
locity dispersion of 0.2 km s -1 ~ 0.2 pc Myr" 1 along the 
baseline of the filament (defined in ij3.2[) . The mean dis- 
tance of sources from the filament was then measured as a 
function of time. The test sources started to show larger 
mean distances than Class I and Flat-Spectrum sources 
after ~ 1 Myr. A mean distance similar to that of Class II 
and Class III objects was reached after ~ 3.5 Myr, which 
is in good agr eement with the evo l ution ary timescale of 
these sources (jLuhman et al.ll2003l [2009h . 

4.2. The Nature of the Dense Cores 
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In W3.21 we denned and examined the dense core popu- 
lation in L 1495-E. To test whether these cores are grav- 
itationally bound en tities, one can est imate the Bonnor- 
Ebert critical mass (|Lada et al.ll2008h 
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Within this framework, in cores with M corc > A/be ther- 
mal motions do not provide enough support against grav- 
ity, and in the absence of other forces they inevitably 
collapse and form protostars. With a mean density of 
riH-) = 1-5 x 10 4 cm~ 3 and a temperature of T = 10 K 
(jTerebev et al.l [2010h . we obtain M BE = 1.5 M Q . Cal- 
culating the Bonnor-Ebert critical mass for all cores al- 
lowed us to analyze the stability of each core separately 
(Fig.[l3]). As a result we find that in L 1495-E, the ma- 
jority of cores is prone to collapse. Without exc eption, al l 
cores that coincide with a dense H 13 CO + core 
exhibit M corc > A/be- Apparently, in these cores the col- 
lapse has already led to a density enhancement, which 
can be traced by H 13 CO + with its critical density of 
n crit ~ 10 5 cm" 3 . 

Compared e.g. to the starless core population in th e 
Pipe Nebula (n H2 = 7.3 x 10 3 cm~ 3 . lLada et al.|[2008h . 
the cores in the Taurus filament are clearly denser by a 
factor of 2. For the cores in the Pipe Nebula, the Bonnor- 
Ebert critical mass is ~ 2M (assuming T = 10 K), 
and only a few cores exceed this value and are prone 
to collapse. Cores in Taurus are also considerably more 
clustered than those in the Pipe Nebula. The nearest 
neighbor separation and two-point correlation function 
(Fig. [7]) show a clear peak at separations of the order of 
0.13 pc, whereas cores in the Pipe Neb ula show separa- 
tions of 0.38 pc ([Rathborne et al.1l2009f) . Fragmentation 
in Taurus happens mainly only on a very distinct spatial 
scale. In the case of thermal fragmentation, such a pref- 
erential length scale would be determined by the Jeans 
length 

/ \ 1/2 

A, = c(^-J , (5) 

where c s is the local sound speed, G the gravitational 
constant, and p the mass density. Assuming c s = 
0.19kms _1 (T = 10 K) and the mean mass density of 
our sample of cores (p = 5.8 x 10 _26 gcm -3 ) we derive 
A ,/ = 0.18 pc, which is only slightly larger than the peak 
in the nearest neighbor distribution of 0.13 pc. This is in 
contrast to the Pipe Ne bula with a separation of 0.38 pc 
(jRathborne et a l. 2009), which exceeds its local Jeans 
length (Aj ~ 0.25 pc). 

We note that due to the unknown inclination angle i we 
can only derive the projected nearest neighbor separation, 
but the true separation between cores remains unclear. 
For i ~ 45° the true nearest neighbor separation would 
match the local Jeans length. VLBA observations al- 
low accurate distance measurement s of selected stars in 
the Taurus Molecular C loud exist (jLoinard 
iTorres et al.ll2007l [2009) , but do not allow to estimate of 
the true spatial orientation of the filament. 

We note, that non-zero inclination angles would not 
only affect the derived nearest neighbor separation, but 
also the aspect ratio. Randomly oriented cores would 
yield the observed mean aspect ratio of R — 1.5, if the 



true aspect ratio is i?t m c ~ 1.7 in the case of prolate, and 
.Rtrue ~ 2.2 for oblate cores, respectively (|Mvers et al.l 
1 19911) . H owever, from a large set of dense cores and Bok 
globules, iRvdenl (|1996f ) argues that cores are generally 
rather prolate than oblate. This shap e is furthermore 
suppor ted by theoretical simulations (jFiege fc Pudritd 
2000b), where helical magnetic fields are able to pro- 
duce a population of prolate dense cores through fila- 
mentary fragmentation. The aspect ratio s of the s e core s 
are in good agreement with the sample of IRvdenl (l996), 
who finds 7? = 2.0. ..2.5. Therefore, under the assump- 
tion of prolate cores, we find cores, which are slightly 
more spherical (i?t rU c ~ 1-7). The same trend is visible, 
when comparing th e Taurus cores from the sample of 
I Lee fc Mversl (|1999f ) to their galactic average. Whether 
this is truly the case, or simply a projectional effect can- 
not be determined. Due to the core formation in the 
filament, we do not expect our sample of cores to be ori- 
ented randomly, but rather inheriting some preferential 
direction. However, the meandering nature of the fila- 
ment does not allow to determine a common inclination 
angle. Considering these limitations, the elongation of 
our cores agrees reasonably well with theoretical predic- 
tions and other observations. 

5. SUMMARY 

In this paper, we study the star formation in the Tau- 
rus filament L 1495. In particular we present a NIR ex- 
tinction map of this filament with unprecedented dynam- 
ical range (Ay ~ 1.5... 35 mag) and resolution (0.9'). 
The main results of this paper are as follows: 

1. The extinction map of L 1495 reveals the highly 
fragmented nature of the filament, harboring a 
population of dense cores preferentially separated 
by the local Jeans length. The mass-per-length in 
the filament is Mij no = 17M pc _1 , indicative of 
its star-forming potential. We find, that the part of 
the filament that harbors no YSOs, namely B 211, 
shows an internal structure different from that of 
the rest of the filament. We argue, that B211 is 
still younger than other parts of the filament, being 
still in process of filament fragmentation. Given its 
high mass reservoir, star formation will inevitably 
ensue. 

2. The dense core population in L 1495-E was inves- 
tigated in detail. A total of 39 dense cores with 
masses between M corc = 0.4 ... 10 M and densi- 
ties of nn 2 = 1-5 x 10 4 cm~ 3 were found. The ma- 
jority of these cores exceeds the critical mass for 
collapse Mbe, and is therefore prone to collapse. 
The high mass tail of the DCMF can be fitted with 
a power-law with an exponent T — 1.2±0.2, a form 
commonly observed also in other star-forming re- 
gions. 
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